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Abstract. We evaluate the transition from zero-sound to first-sound behaviour with 
increasing collisionality in the propagation of density waves through an ultracold 
gaseous mixture of fermionic atoms confined in the normal state inside a cigar-shaped 
harmonic trap. We study for this purpose the evolution of the one-body distribution 
functions associated with a density perturbation generated in the central region of the 
cloud, as obtained by solving numerically the Vlasov-Landau equations. We examine a 
variety of trap anisotropies and of repulsive or attractive interaction strengths between 
the components of the mixture, and the speed of propagation of the density disturbance 
is found to decrease in both cases as the magnitude of the coupling strength is increased. 
The results are compared with the values of the speed of zero sound and of first sound, 
as obtained analytically from the limit of vanishing collisionality and from linearized 
hydrodynamics. The main effects of the quasi-one-dimensional confinement are the 
stabilization of zero-sound excitations in the attractive regime before collapse and 
the lowering of the hydrodynamic sound velocity by a factor y/3/5 relative to three- 
dimensional behaviour. 
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1. Introduction 

Several experiments on mixed Fermi gases consisting of two components with different 
pseudospins have probed their collisional properties both in the normal [f] and in 
the superfluid state [2, 3]. In the experiments performed by the JILA group, the 
collisionality of a 40 K mixture in the normal state was tuned by varying either the 
atomic density [1,4] or the off- resonant value of the scattering length [5,6]. Gensemer 
and Jin [1] examined the frequency and damping of the dipolar modes of the mixture 
as functions of its collisionality, showing that, while in the collisionless limit the two 
components oscillate independently with vanishing damping, on increasing the strength 
of the mutual repulsive interactions the damping first increases and then decreases 
until in the hydrodynamic limit the two components oscillate together at the same 
frequency without damping. The same behaviour has been obtained from numerical 
simulations [7,8]. The absence of damping in the two opposite limits signals the existence 
of long-lived collective sound-wave excitations in suitable system geometries, the so- 
called zero sound in the collisionless limit and the hydrodynamic first sound in the 
strongly collisional regime. 

More recently attention has focused on sound-wave propagation in a Fermi gas with 
attractive interactions ranging in strength from the BCS to the BEC regime [9-12]. 
Below the critical temperature for the transition to a superfluid state, the presence 
of both normal and superfluid components leads to two acoustic modes corresponding 
to the propagation of density waves and thermal waves [9, 13]. These modes contain 
information on the equation of state of the gas and on the geometry of the system [1 1] 
and impose constraints on the appropriate theoretical models to describe the BCS-BEC 
crossover [10]. Above the critical temperature the speed of propagation of a density 
perturbation still gives information about the shape of the equilibrium density profiles 
and about the collisionality of the gas relative to the frequency of the perturbation. 

The behaviour that we have recalled above for the attenuation of sound waves 
as a function of collisionality was predicted in the early work of Landau on normal 
liquid 3 He [14]. In the experiments of Abel et al. [15] on this liquid, the collisionality 
was explored by tuning the frequency of an acoustic disturbance. The experiment 
demonstrated that the transition from collisional to collisionless dynamics is signalled 
by a drastic change in the frequency dependence of the attenuation and showed that 
the ratio of the speeds of zero and first sound is very much smaller than the value v^3 
expected for a homogeneous fluid at weak coupling. In a dilute Fermi gas one should 
expect that the interactions introduce instead only slight modifications to the bare sound 
velocities in the ideal Fermi gas. 

The main purpose of the present work is to provide direct theoretical evidence 
for the transition from zero to first sound in a normal Fermi gas by evaluating the 
propagation of density waves in a two-component mixture confined inside an elongated 
cigar-shaped harmonic trap. We have in mind an experiment in which the collisionality 
could be varied by tuning the scattering length a as can be realized by exploiting 
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Feshbach resonances [5] and the density perturbation could be generated by a laser 
beam focused at the centre of the trap [16]. We perform numerical simulations of the 
dynamics of the gas, in which collisions are taken into account by solving the Vlasov- 
Landau equations (VLE) for the one-body distribution functions of the two fermionic 
components. We observe that on increasing the scattering cross-section the velocity 
of propagation of the density distortions, which is given by the Fermi velocity vp in 
the collisionless limit, diminishes for both attractive and repulsive coupling. If the 
mean-field corrections are negligible, the hydrodynamic limit for the confined gas in the 
elongated trap yields a propagation velocity v-p/yh instead of Vf/V3 as appropriate to 
a homogeneous gas. Most remarkably, the zero-sound excitation is found to be stable 
even at weak attractive couplings, before the attractions drive the gas to collapse. 

The paper is organized as follows. In Sec.EJwe introduce the fermion mixture under 
study and discuss the dynamical regimes that can be realized in a dilute gaseous state, 
a detailed analytic calculation for zero sound in a cylindrical trap being reported in an 
Appendix. In Sec. EJwe present our treatment of a mixture with finite collisionality and 
discuss our numerical results for the velocity of propagation of density waves. Finally, 
Sec. |U offers a summary and the main conclusions of the work. 



2. The model system 

We consider a mixture of two species of fermionic atoms of mass m confined in a 
harmonic trap of the form 

Kxt(r) = im^i(ri + AV), (1) 

where A = uo z /uo± is the anisotropy parameter. The atoms of each species are spin- 
polarized and the gas is highly dilute, so that at very low temperatures only interactions 
between atoms of different species are allowed. The gas is described in terms of the one- 
body distribution functions f^(r,p,t) for each species s and the particle densities are 
obtained by performing an integration over momentum space, 

n W (r,*) = / (^/ W ( p >P»*)- ( 2 ) 
The distribution functions at equilibrium are given by the local Fermi-Dirac distributions 



fo\r,p) = < exp 



(3) 



where (3 = 1/ksT with T the temperature, //W is the chemical potential ensuring the 
normalization N s = J f^ s \r,p) d?r d 3 p/h 3 , and the mean-field potentials read 

U&(r) = V ext (r)+gn®(r). (4) 

Here s denotes the component different from s and g = AixTi 2 a/m is the coupling strength 
between atoms of different species, with a being their s-wave scattering length. 

The calculation of the equilibrium density profiles no requires the self-consistent 
solution of equations (jSJ)-® for each set of system parameters. At T = the results 
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(s) 

of the Thomas- Fermi approximation are recovered, giving n as the solution of the 
equation fj,(n^) + V cxt (r) = /j,( s \ with fi(n) the equation of state of the gas 

^ ) ) = ^{^ 2 n^f" + gnf. (5) 

We shall study a mixture of 2000 40 K atoms at T = 0.2 Tp, Tp being the Fermi 
temperature. The atoms are equally distributed in two spin states and confined inside 
a trap with radial frequency u± = 100 s _1 . Hereafter we will remove the indices s and 
s when there is no ambiguity and set /i^ = ^ = \i and = = n. 



2.1. Sound propagation 

The time evolution of the one-body distributions strongly depends on the frequency 
of the inter-species collisions. When this is much lower than the trap frequency, 
the evolution of the gas is governed by deformations of the Fermi sphere and the 
dynamics is termed collisionless. This is the regime of zero sound. In a homogeneous 
mixture the collisionless sound velocity is given by Co = r/vp, with the Fermi velocity 
vp = [2(/x — gn^/m} 1 / 2 . The ratio rj can be calculated in Landau's Fermi-liquid theory 
describing the mean-field interactions through appropriate Landau parameters. For a 
two-component Fermi gas with negligible interactions between atoms of the same spin, 
rj satisfies the relation [14, 17] 

^2 + 1-1 = -^-. (6) 
2 i] — I 2mvp a 

In a trap one may use a local-density approximation and assign the position-dependent 
Fermi velocity vp(r) = [2(/i— V cxt (r)— gno(r)) /m] 1 / 2 to the atoms, provided the variation 
of the density profile is sufficiently smooth. In this way the confinement is accounted 
for only through the position dependence of the density profiles. 

However, a very tight confinement may drastically modify the density of states of 
the gas and thus the dynamic of density perturbations. The zero sound velocity c 
in strongly elongated traps (A — > and V cxt (r) — > muj 2 _rj_/2), in a regime where the 
degrees of freedom in the azimuthal plane are frozen by the confinement, is given by 
(see Appendix) 

o {, V /2 
co = vUl + — o , 7 

with Vp being the local Fermi velocity at the centre of the trap and 2ira 2 L the radial 
section of the atomic cloud. This velocity coincides with that of a strictly one- 
dimensional (ID) Fermi gas [18,19] with an effective coupling 2h 2 a/(ma 2 L ) embodying 
the transverse confinement [20]. 

Interactions between the two components make the sound velocity deviate from 
Vp by only a few percent for couplings up to a = 2 x 10 4 Bohr radii for both the 
homogeneous and the cigar-shaped gas [17]. The effects of the interactions have also 
been investigated in the random-phase approximation for a Fermi gas in a spherical trap 
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in the collisionless regime [21,22]. The main outcome is the appearance of a fragmented 
excitation spectrum at strong couplings and of Landau-damped modes for attractive 
interactions. However, the confinement introduces a discrete single-particle excitation 
spectrum that may weaken the strength of the damping. In a cigar-shaped confinement, 
the sound velocity given in equation (JZJ) is well defined for repulsive interactions as well 
as sufficiently weak attractive interactions, analogously to zero sound in a ID Fermi 
gas [19]. 

When the dynamics is instead governed by the hydrodynamic equations, the (first) 
sound velocity in a homogeneous gas is given by 

v F ( 2amv F \ 1/2 

showing that the explicit dependence on the scattering length is qualitatively similar 
to that in equation (JJJ). In addition to the density modes, an oscillation of the spin 
density £ = n^ s ' — may develop depending on how the gas is perturbed [23-25] . In 
fact, for strong repulsions the two species tend to segregate into two non-overlapping 
regions. For attractive interactions, the system will become unstable and experience a 
sudden collapse when the sound velocity vanishes. However, at very low temperatures 
the collapse may be prevented by a transition to a superfluid phase [2, 26, 27] or by 
pairing of the atoms into molecular dimers [28-30]. 

With the purpose of examining the confinement effects on the propagation of first 
sound, we consider again a radially confined Fermi gas and assume a density wave 
travelling in the axial z direction. In this case we have shown [11] that the first sound 
velocity C\ depends on the equation of state through the relation 

ci = Q J n d 2 r ± / J {dv/dn\ no y l d 2 r^j 12 (9) 

where the integration is on the azimuthal plane. Equation (jHJ) reduces to equation (JBJ) 
for a homogeneous Fermi gas and yields C\ = v F / \/h for the cigar-shaped Fermi gas 
without mean-field interactions. More generally, numerical integration can be carried 
out once the equilibrium density is known. 

3. Density waves at finite collisionality 

Neither the collisionless nor the hydrodynamic approach is appropriate when the number 
of collisions is sizable but not sufficient to establish local equilibrium in the gas. One 
must then resort to kinetic equations and can use a Boltzmann equation method to take 
into account binary collisions in a semiclassical manner. The evolution of the one-body 
distributions is determined by the VLE 

d t f {s) + - • V r /W - V r E/« • V p /W = C^if( s \ /«], (10) 
m 

where the collision integral C^ s ' reads 

cW = s^sf /*» ia ' v 1(1 " /W)(1 " 
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-/ (S) /2 (5) (l-/3 (S) )(l-/f)]. (11) 

Here = f( s \r,p,t) and ff' = f( s \r,Pi,t), dVtj is the element of solid angle for 
the outgoing relative momentum p 3 — p 4 , o = Arc a 2 is the scattering cross- sect ion, and 
v = \v — v 2 \ is the relative velocity of the incoming particles. The collision satisfies 
conservation of momentum (p + p 2 = P3 + Pa) and energy (e + e 2 = £3 + £4), with 
Si = p 2 /2m + U^\ 

The propagation of a density perturbation is evaluated by performing the following 
numerical experiment. First we construct the equilibrium density profiles according 
to equations (J2J)-® for a given confinement configuration and interaction strength. 
A Gaussian potential U exp(— z 2 /2w 2 ) with amplitude U = 10hu± and width w = 
(U/ (muj^)) 1 / 2 is then switched on at the centre of the trap to simulate a laser beam 
expelling particles from the centre of the cloud. At later times these density distortions 
are seen to move towards the ends of the trap. Since each perturbation has velocity 
mainly in the z direction, we focus on the evolution of the integrated density profile 
n(z) = J n(r)d 2 r±. Experimentally, this density can be obtained directly from in-situ 
measurements of the three-dimensional density profiles. In figure Q we plot n(z) at 
several times for a trap with A = 0.1 and for three values of the scattering length, 
namely a = —1.5 x 10 4 , 80, and 2 x 10 4 Bohr radii ao- At larger attractions the gas 
collapses, whereas for larger repulsions the effects of demixing become noticeable. We 
follow the positions of the main density peaks created by the Gaussian potential and 
plot in figure El their evolution in time. After a hole in the density has been created, we 
observe that density pulses travel with constant speed in the neighbourhood of the trap 
centre. For a = —1.5 x 10 4 ao and 2 x 10 4 ao it can be seen that the slope of the peak 
positions clearly departs from cq and moves closer to c±, as the hydrodynamic regime is 
being approached. 

The results for the sound velocity as a function of the coupling are summarized in 
figure EJ where we show the mean value (v z ) = ^({v^ s ') + (i^ )) of the time-averaged 
velocities of the two components of the mixture as obtained from the solution of the 
VLE. These have been evaluated over a time interval At = 1.5/ u± starting at t = l/^_i_- 
In the same figure we also report the results of the analytical models introduced in 
Sec. 12 .1\ i.e. the zero-sound velocity c in a strongly elongated trap (cf. equation (J7J)), 
and the first-sound velocity c\ in a strongly elongated trap as given in equation © 
with = n^ s \ In the latter case we have used the zero-temperature equation of state 
(0) and set fi to the value obtained in the VLE simulation. By solving the VLE in 
the absence of collisions we also find that the speed of the density perturbations is well 
defined for both positive and negative couplings and follows the trend of c as a function 
of a (not shown in the figure). We emphasis that in the absence of collisions zero sound is 
a well defined collective excitation in the gas under quasi-lD confinement even when the 
interactions between the particles are attractive, contrary to the case of a homogeneous 
3D gas where these excitations are stable only for interparticle repulsions [14]. On the 
other hand, in the case of repulsive interactions our results for cq are very close to those 
obtained from equation for the homogeneous gas. 
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Figure 1. Time evolution of the integrated density profile n(z) as a function of axial 
position z (in units of a± = y/h/(mu)±)) for 20 40 K atoms in a trap with A = 0.1. 
The left, middle, and right panels correspond to a — —1.5 x 10 4 , 80, and 2 x 10 4 Bohr 
radii, respectively. The bottom profile refers to the initial equilibrium state and the 
subsequent ones are separated by time intervals of 0.5/ui± = 5 ms. The later profiles 
have been vertically displaced for the sake of visibility. 



From the full solution of the VLE including collisions, we observe that at low 
coupling the velocity of the density perturbations is close to cq within the error bars. In 
this regime the value of the sound velocity is dominated by the mean-field interactions 
and their effect on sound propagation can be understood as mainly resulting from 
changes in the local density: in particular, for attractive coupling the Fermi velocity 
and the sound velocity rise as the density of the gas is increased. At stronger coupling 
strength the sound velocity drops for both positive and negative coupling, as opposed to 
the behaviour of the collisionless velocity Cq. This drop is more pronounced for attractive 
coupling as the frequency of collisions increases with \a\ as well as with the gas density, 
whereas at strongly repulsive interactions the collisions increase more slowly as the two 
components approach spatial demixing. This effect is not included in the calculation of 
c\. At strongly attractive interactions the frequency of collisions is much larger than 
u z and therefore the gas is strongly collisional. The difference between (v z ) and c\ 
at a = —1.5 x 10 4 ao may be attributed to the appreciable amplitude of the density 
perturbations as well as to the confinement in the axial direction. At strong coupling, 
although the rate of collisions in the axially confined inhomogeneous cloud varies at 
each point in space during the propagation of the pulse, the sound wave travels without 
observable damping in the neighbourhood of the trap centre. 

The behaviour of (v z ) for two values of the anisotropy parameter at fixed u± is 
illustrated in figure 0] In the top panel we report the Fermi velocity at the cloud centre 
for both components, showing that in the more elongated cloud the density of the gas 
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Figure 2. Location z p /a± of the density peaks in n(z) as a function of time (in units 
of 1/oj±). The top, middle, and bottom panels correspond to a = —1.5 x 10 4 , 80, and 
2 x 10 4 Bohr radii. The dashed and dotted lines refer to propagation at the zero-sound 
velocity Co and at the first-sound velocity c\, respectively. 
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Figure 3. Average velocity (v z ) (in units of y / Tiuj±/m) as a function of the scattering 
length a (in units of the Bohr radius ao) for A = 0.1. The thick solid line and the 
dashed line show the sound velocity calculated in the absence of damping as for first 
sound and zero sound, respectively. The dots with error bars are the results of the 
VLE simulation, the errors being estimated from the noise in the slope of the density 
peak positions in figure |3 The thin solid line is just a guide to the eye. The scale at 
the top gives the quantum collision rate T q (in units of uj z ). 



and thus the Fermi velocity are lower. Therefore, at a fixed value of a the collisionality 
of the gas with A = 0.06 is lower and (v z )/vp is higher (bottom panel of figure 0J. On 
the other hand, if we fix the average trap frequency and thus the density at the centre, a 
deformation of the trap would increase the number of collisions occurring during sound 
propagation and therefore facilitate the transition towards the hydrodynamic limit [31]. 
The effect of demixing at strong repulsive interaction is also illustrated in the top panel 
of figure |U 

4. Summary and concluding remarks 

We have studied the propagation of sound waves in a mixture of fermionic atoms confined 
inside a strongly elongated cigar-shaped trap and exposed to the role of collisions in 
determining the speed of propagation. We have compared our results from the numerical 
solution of the Vlasov-Landau equations with two limiting behaviours: zero sound and 
first sound propagating in an infinitely long cylindrical cloud with a very tight transverse 
confinement. We have found substantial agreement between the velocity of density 



Sound wave propagation in strongly elongated fermion clouds at finite collisionality 10 

5 I 1 1 1 1 1 1 — 

A = 0.1 — •— 

A = 0.06 




_i i i i i i 

10 4 10 4 2 x 10 4 3 x 10 4 4 x 10 4 

a/a 




i i i i i i 

-10 4 10 4 2 x 10 4 3 x 10 4 4 x 10 4 

a/a 



Figure 4. Sound velocity in elongated traps with A = 0.1 or A = 0.06. Top panel: 
Local Fermi velocity Wp at the cloud centre (in units of yjTiuj^/m) as a function of a 
(in units of the Bohr radius). Bottom panel: Average sound velocity (in units of u F ) 
as a function of a (in units of the Bohr radius). 
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perturbations in the collisionless regime and in the highly collisional regime with those 
of zero sound and first sound, respectively. The differences may be attributed to the 
finite amplitude of the density perturbations that we have produced and to the axial 
confinement. 

It would be interesting to study the effect of the presence of a superfluid component 
on the zero-to-first sound transition at finite temperature. This can be done by using 
suitable transport equations for a trapped Fermi gas below the superfluid transition 
temperature, as already treated by Urban and Schuck [12]. 
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Appendix A. Zero sound in a cylindrical trap 

We apply the transport equation (fTUj) with = and V ext = muj 2 1 r\/2 = V to 
evaluate the propagation of zero sound in a tight cylindrical trap. A density perturbation 
travels mainly along the axis of the cylinder when the radial confinement is strong, and 
the distribution function can be written as 

f is) (r,p) = ti s \r ± ,p) + ft\r ± ,p)e^-^ (A.l) 
where /q = 6(p^ — J~ — V(r±) — gnQ\r±)) is the equilibrium distribution. Around 
the centre of the azimuthal plane one can write = 9{pS 8 ^ — Jm — fm — \rn£j 2 _r 2 1 ) 
with #W = fi^ - gn { o\o) and u± = u±.(l - Bgnjp (0) /p/i)) 1 / 2 . Only the first level 
in the radial harmonic potential with frequency ti± is occupied by the fermions in 
a very tight radial confinement, and the distribution function takes the simple form 

/o (s) = 0(£ (s) -fi-^). 

The terms in the transport equation involving the derivatives with respect to rj_ 
and p± disappear if f[ s ' is a function of the total energy analogously to . Equation 
(|TU]l then becomes 

ftp. \ fi , ) _k_P 1 a_g> [ jp i , ) = ^ 



\ m J m de J (27rh) 3 ' 

where we have set e = pi/ (2m). The form of equation (jA.2|) suggests to write f± as 
proportional to udf^ /de, v being defined as 

-— if kp z > 

—— if kp z < 
l + i] 
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with rj = u/(kv^ ) and = \/2(/z( s ) — hQ±)/m. Thus equation (jA.2|) can be rewritten 



as 



1 + r^/(^' !, -f|- K ^ - < A - 4 » 



Finally, we obtain 
V = 



\ 



1 + — ; M ■ ( A -5) 



~2 (5) 

ira , mvi? 



by integrating over all p z and p± G [0, y/2hmu±]. This is equation ((7j) in the main text. 
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